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ABSTRACT 

We have observed the S(0), S(l), S(2), S(3), S(4) and S(5) rotational hnes 
of molecular hydrogen (H2) towards the peak of the photodissociation region 
(PDR) associated with the reflection nebula NGC 7023. The observed H2 hne 
ratios show that they arise in warm gas with kinetic temperatures ~ 300 - 700 
K. However, the data cannot be fitted by an ortho- to para- (OTP) ratio of 3. 
An OTP ratio in the range ~ 1.5 - 2 is necessary to explain our observations. 
This is the first detection of a non-equilibrium OTP ratio measured from the H2 
pure-rotational lines in a PDR. The existence of a dynamical PDR is discussed 
as the most likely explanation for this low OTP ratio. 

Subject headings: ISM: abundances — ISM: individual (NGC 7023) — reflection 
nebulae — stars: individual (HD 200775) — stars: pre-main-sequence — 
infrared: ISM: hnes and bands 
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1. Introduction 

The Short Wavelength Spectrometer (SWS) on board the Infrared Space Observatory 
(ISO) has allowed us, for the first time, to observe the pure rotational spectrum of H2 and 
study the OTP ratio in warm (T~ 100 - 1000 K) regions. Previous to ISO observations, 
the OTP ratio of H2 had been estimated towards regions heated by UV radiation and 
shocks to temperatures > 2000 K using the relative strengths of the vibrational lines. In 
outflow regions where shock excitation is the main heating mechanism, observations of the 
vibrational lines typically reveal OTP ratios of ~ 3 (Smith, Davis & Lioure 1997). This 
is the expected value in thermodynamic equilibrium for gas with temperatures > 200 K. 
However, towards the regions heated mainly by the UV radiation (PDRs), the vibrational 
lines give OTP values in the range 1.5 - 2 (see e.g. Chrysostomou et al. 1993). The pure 
rotational hues of H2 have been detected by ISO toward several Galactic regions (S140: 
Timmermann et al. 1996; Cepheus A West: Wright et al. 1996; BD+40°4124: Wesselius 
et al. 1996; HH54: Neufeld et al. 1998). The observations towards all these regions 
except HH54 are consistent with the H2 rotational lines arising in gas with temperature 
> 200 K and with an OTP ratio of 3. This value is not consistent with those derived 
from the vibrational lines in dense PDRs but is in agreement with theoretical predictions. 
Recently Sternberg & Neufeld (1999) have argued that the low OTP ratios measured 
from vibrational hues do not represent the actual OTP ratio in PDRs but it is simply a 
consequence of the optical depth effects in the fluorescent pumping of the vibrational lines. 
The non-equilibrium OTP ratio measured towards HH54 has been interpreted as arising in 
hot-shocked excited gas that has not reached the equilibrium yet. In this interpretation, 
the observed OTP ratio would be the legacy of an earher stage in the thermal history of the 
gas when the gas temperature was 90 K. 

In this Letter we report the observations of the H2 rotational lines towards the 
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prototypical PDR associated with the reflection nebula NGC 7023. An OTP ratio in the 
range 1.5 - 2 is derived from our observations. This is the flrst detection of an OTP ratio < 
3 in a PDR based on the pure H2 rotational lines. 

2. Observations 

The observations of the S(0), S(l), S(2), S(3), S(4) and S(5) rotational lines of H2 
towards the peak of the PDR associated with NGC 7023 (R.A. (2000): 2l'^01'"32^5, 
Dec(2000): 68°10'27".5). were made using the Short- Wavelength Spectrometer (SWS) (de 
Grauuw et al. 1996) on board the Infrared Space Observatory (ISO) (Kessler et al. 1996) 
during revolution 514 with a total on-target time of 5149 s. At the spectral resolution of 
this mode (A/AA ~ 1000-2000) all the observed lines are unresolved. Data reduction are 
carried out with version 7.0 of the Off Line Processing routines and the SWS Interactive 
Analysis at the ISO Spectrometer Data Center at MPE. Further analysis has been made 
using the ISAP software package. The uncertainities in the calibration are of 15 %, 25 %, 
25 %, 25 %, 20% and 30% for the S(5), S(4), S(3), S(2), S(l) and S(0) lines respectively 
(Salama et al. 1997). 

3. H2 rotational lines 

In Fig. 1 we present the spectra of the S(0), S(l), S(2), S(3), S(4) and S(5) H2 lines 
towards the peak of the PDR associated with NGC 7023. The observed intensitites and 
some interesting observational parameters are presented in Table 1. The data have been 
corrected for dust attenuation using the extinction curve of Draine & Lee (1984), and a 
value of 0.43 for the dust opacity at 0.55 fj,m. This value has been derived from the ISO 
LWSOl spectrum (Fuente et al. 1999, hereafter paper II). The extinction for the S(0), 
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S(l), S(2), S(3), S(4) and S(5) lines is very small and amounts to 0.01, 0.02, 0.02, 0.04, 
0.02 and 0.01 respectively. Fig. 2 shows the rotational diagram for the H2 corrected for 
extinction effects and assuming extended emission. (Note that the errors in Fig. 2 are 
entirely dominated by the calibration uncertainities.) 

This diagram shows that the ortho-H2 levels have systematically lower N„/gu values 
(where and g^ are the column densities and degeneracies of the upper levels of 
the transitions) than the adjacent J-1 and J+1 para-H2 levels producing a "zig-zag" 
distribution. In fact, the ortho-H2 levels seem to define a curve which is offset from that of 
the para-H2 levels (see Fig. 2). This is the expected trend if the OTP ratio is lower than 3. 
The offset between the two set of data is systematic and seems to show the trend of being 
larger for the low energy levels (S(0) and S(l) transitions), than for the high energy levels 
(S(2), S(3), S(4) and S(5) transitions). 

The offset between the ortho- and para- H2 curves is larger than the observational 
errors and cannot be due to the different apertures for the different lines. For the case of a 
point-like source we would have to correct by a factor of 1.93 the intensity of the S(0) fine 
and by a factor of 1.35 the intensities of the S(l) and S(2) fines. In this case, the offset 
between the para- and ortho- curves would increase. Based on the spatial distribution of 
the HI(21 cm) fine (Fuente et al. 1998) and those of the CII (158 //m) and 01 (63 /xm) 
lines (Chokshi et al. 1988, paper II), we will assume in the following a beam filfing factor of 
1 for all the observed H2 lines. Although the exact shape of the rotational diagram depends 
on the excitation conditions of the region, the "zig-zag" features cannot be explained 
by any model which assumes an equilibrium OTP ratio. In thermodynamic equilibrium, 
the rotational temperature between an ortho- (para-) H2 level and the next J-l-1 level 
should always increase with the energy of the upper level giving rise to a smooth curve 
in the rotational diagram regardless of the temperature profile of the region. A "zig-zag" 
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distribution implies that for some pair of levels the rotational temperature decreases with 
the energy of the upper level. One can only get this effect with a non-equilibrium OTP 
ratio. 

We have compared our data with the models by Burton et al. (1992) in order to 
estimate the incident UV field, density and OTP ratio. Burton et al. (1992) assumed a 
constant value of 3 for the OTP ratio in their calculations. The best fit to our data is for 
an incident UV field of Go = 10^ in units of the Habing field and a density of n = 10^ 
cm~^, but it systematically underestimates the intensity of the para-H2 transitions and 
overestimates the intensities of the ortho-H2 transitions (see open triangles in Fig. 3). 
This is the expected behavior if one assumes an OTP ratio of 3 and the actual OTP ratio 
in the region is < 3. We have corrected the line intensities predicted by the model for 
different values of the OTP ratio assuming that the total amount of H2 molecules at a given 
temperature and the line ratios between levels of the same symmetry are not affected by 
the OTP ratio. These assumptions are correct for the low- J transitions where the coUisional 
excitation dominates. In Fig. 3 we compare our data with the predicted diagram for an 
OTP ratio of 3 which correspond to the model without any correction (open triangles), 2 
(open circles) and 1 (open squares) . The diagram for an OTP ratio of 3 clearly does not fit 
any of our observational points. We find the best fit to our data with an OTP ratio of 2. In 
this case, the predicted intensities for the S(l), S(2) and S(3) lines are in agreement with 
the observed values. Only the predicted intensity for the S(0) line is not consistent with the 
observations. To fit the intensity of the S(0) line it is necessary to assume an OTP ratio of 
< 1.5 (see the case OTP = 1 in Fig. 3), but in this case, we will have a worse fit for the 
S(l), S(2) and S(3) hues than for an OTP ratio of 2. As discussed in Section 4, this suggests 
a possible variation of the OTP ratio with the temperature. The cooler gas emitting only 
in the S(0) line could have an OTP ratio of < 1.5, while the warmer gas emitting in the 
S(l), S(2) and S(3) lines could have an OTP ratio of 2. However, to establish this variation 
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unambigously, it is necessary to have a very accurate knowledge of the excitation conditions 
in the region. To be conservative, we will conclude that the OTP ratio is in the range 1.5 - 

2 in this PDR. Martini et al. (1997) found an OTP ratio of 2.5±0.4 using the near-IR H2 
vibrational lines towards their position 1 which is 20" offset from ours. Taking into account 
that because of the optical depth effects in the vibrational lines this value is just a lower 
limit to the actual OTP ratio, it proves that the OTP ratio is close to 3 for the gas with 
kinetic temperatures T^ > 2000 K. This difference between the OTP ratios derived from 
the rotational and vibrational lines argues in favor of a variaton of the OTP ratio with the 
kinetic temperature. 

Since the OTP ratio in this source is different from the equilibrium value, the rotation 
temperature between an ortho- and a para- level does not represent an estimate of the 
gas kinetic temperature. To estimate the gas kinetic temperature we have calculated the 
rotation temperature between levels of the same symmetry. For the para-H2 levels, we have 
derived a rotation temperature of ~ 290 K from the S(0) and S(2) lines and of ~ 500 K 
from the S(2) and S(4) lines. For the ortho-H2 levels, the derived temperature is ~ 440 from 
S(l) and S(3) lines and ~ 700 K from the S(3) and S(5) lines. Based on these calculations 
we conclude that the OTP ratio is ~ 1.5 - 2 in the gas with kinetic temperatures ~ 300 - 
700 K. We derive a total H2 column density of 5 10^° cm~^ assuming a rotation temperature 
of 290 K between the J = and J = 2 para-levels and of 440 K between the J = 1 and J = 

3 ortho-levels. 

4. Discussion 

Based on the H2 rotational hues data we have derived an OTP ratio in the range of 
1.5 - 2 for the gas with kinetic temperatures ~ 300 - 700 K towards the reflection nebula 
NGC 7023. This is the second object with a non-equilibrium OTP ratio measured from the 
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H2 pure rotational lines. The first non-equilibrium OTP ratio was detected towards the 
outflow source HH54, and interpreted as arising in shock-heated gas that has not reached 
the equilibrium yet. This interpretation is not plausible for NGC 7023. The high dust 
temperature,Td~ 40 K, and the detection of the Sill (34.8 /xm) and CII (157.7 //m) hues 
(paper II) proves the existence of a PDR in this region. Although a shock component might 
also exist (sec Martini et al. 1997), the heating is dominated by UV photons. It is then 
expected that most of the warm H2 arises in the PDR. 

A non-equilibrium value of the OTP ratio is not expected for the physical conditions 
prevailing in a dense PDR. The initial OTP ratio after the H2 formation is very uncertain. 
Because of the large exothermicity of the formation process, it is expected to be 3 . 
However, this OTP ratio can change if the H2 molecules remain on the grain surface long 
enough to reach the equilibrium at the dust temperature. In this case, the OTP ratio after 
H2 formation will be that of the equilibrium at the grain temperature. After the ejection 
of the H2 molecules to the gas phase, exchange reactions with H and H"*" change the OTP 
ratio until achieving the equilibrium at the gas temperature. For the gas temperatures 
traced by the H2 rotational lines (> 200 K), the equilibrium OTP ratio is 3. One possibility 
to explain the low OTP ratio observed in NGC 7023 is to suppose that the OTP ratio in 
the H2 formation is lower than 3, and once in the gas phase, the H2 molecules are destroyed 
before attaining the equilibrium value at the gas temperature. The dust temperature in 
a PDR with Go = 10^ and n = 10^ cm~^ is < 100 K. In particular, the dust temperature 
measured towards the PDR peak in NGC 7023 using our LWSOl spectrum is 40 K which 
corresponds to an equilibrium OTP ratio of ~ 0.1. The OTP conversion in the atomic 
region is dominated by H - H2 collisions with a rate of 10~^^n s~^ where n is the atomic 
hydrogen density (see Sternberg & Neufeld 1999 and references therein). The hydrogen 
density derived from the HI image published by Fuente et al. (1998) is ~ 5 10^ cm~^. A 
similar value (~ 0.5 - 1 10^ cm~^) was derived by Chokshi et al. (1988) based on the 01 and 
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CII lines. However, densities larger than 10^ cm~^ are derived from molecular data (Fuente 
et al. 1996, Lemaire et al. 1996, Gerin et al. 1998). The density of atomic hydrogen in the 
region with ~ 300 K is expected to be about an order of magnitude lower than the H2 
density. Then we assume an atomic hydrogen density of ~ 10^ cm~^ in our calculations. 
With this value, the OTP conversion rate due to H2 - H coUisions is 10~^ s~^. At the cloud 
surface the unshielded photodissociation rate is ~ 5 lO^^^Go ~ 5 10^^ s~^, i.e., 2 orders 
of magnitude larger than the OTP conversion rate. Then, in the cloud surface, before 
self-shielding is important (A^ < 0.3 mag, Tfe> 500 K), an OTP ratio lower than 3 can be 
explained by assuming that the OTP ratio in the H2 formation is the equilibrium value at 
the grain temperature. Deeper into the molecular cloud, when H2 is self-shielded, one can 
consider that the H2 destruction rate is similar to the H2 reformation rate which is given 
by ~ 3 10"^^n s~^. In this region, a lower limit to the OTP conversion rate is given by the 
conversion rate due to H2 - collisions which is ~ 10^^ ''n s^^ (Sternberg & Neufeld 1999). 
The OTP conversion rate is of the same order than the destruction rate and the OTP ratio 
is expected to be close to the equilibrium. According to these estimates one expects to have 
an OTP ratio close to 3 at low temperatures and a non-equilibrium OTP ratio at higher 
temperatures. The contrary trend is derived from our observations. The high energy S(l), 
S(2) and S(3) hues are fitted with an an OTP ratio of 2 while the low energy S(0) hne is 
better fitted with an OTP ratio of 1.5. 

Another possibihty to explain the non-equilibrium OTP value in NGC 7023 is to 
consider the case of a dynamic PDR, i.e. the dissociation front is advancing into the 
molecular cloud. In this case we do not have to assume a non-equilibrium OTP ratio in the 
H2 formation. The PDR is being fed continuously by the cool gas of the molecular cloud 
in which the equihbrium value of the OTP ratio is lower than 3. This gas is heated by the 
stellar UV radiation to temperatures > 200 K but leaves the PDR before attaining the 
equilibrium OTP ratio at this temperature. In this scenario, the gas which is expected to 
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have an OTP ratio smaller than the equilibrium value is the gas that has more recently been 
incorporated into the PDR, which is also the gas at lowest temperature. The OTP ratio 
is expected to increase and reach values close to 3 at high temperatures. This behavior is 
consistent with the trend observed in our data. As discussed above, the OTP H2 conversion 
is mainly due to H - H2 collisions with a conversion rate of ~ 10~^ s~^. To have a significant 
fraction of the gas with a non-equilibrium OTP ratio the photodissociation front must 
advance about 1-2 mag in the conversion time, i.e., it must penetrate into the molecular 
cloud at a velocity of ~ 10^n~^ kms~^. Skinner et al. (1993) proposed the existence of an 
anisotropic ionized wind associated with this star based on radio continuum observations. 
Later, Fuente et al. (1998) detected an HI outflow with a velocity of ~ 7.5 kms~^. They 
proposed that this outflow is formed when the gas that has been photodissociated by the 
UV radiation, is accelerated by the stellar winds along the walls of a biconical cavity. The 
velocity of the HI outflow cannot account for the velocity at which the photodissociation 
front must advance into the molecular cloud to have an OTP ratio lower than 3 unless the 
OTP conversion rates are severely overestimated (by a factor of > 10) and/or the density of 
the outflow is > 10^ cm~^. In spite of this problem, a photodissociation front penetrating 
into the cloud because of the outflow seems to be the most plausible explanation for the 
low OTP ratio measured in this region. The existence of an outflow can also explain the 
difference between this region and other PDRs hke S140 in which an OTP ratio of 3 has 
been derived from the H2 rotational lines. 

We would like to thank the referee for his/her helpful comments. This work has been 
partially supported by the Spanish DOES under grant number PB93-0048 and the Spanish 
PNIE under grant number ESP97-1490-E. N.J.R-F acknowledges Conserjeria de Educacion 
y Cultura de la Comunidad de Madrid for a pre-doctoral fellowship. 



- 11 - 

REFERENCES 

Burton M.G., HoUenbach D.J., Tielens A.G.G.M., 1992, ApJ 399, 563 
Chokshi A., Tielens A.G.G.M., Werner M.W., Castelaz M.W., 1998, ApJ 334, 803 
Chrysostomou A., Brand P.W.J.L., Burton M.D., Moorhouse A., 1993, MNRAS 265, 329 
Draine B.T. & Lee H.M., 1984 ApJ 285, 89 

Fuente A., Martm-Pintado J., Neri R., Rogers C., Moriarty-Schieven G., 1996, A&A 310, 

286 

Fuente A., Martm-Pintado J., Rodriguez-Franco A., Moriarty-Schieven G., 1998, A&A 339, 
575 

Fuente A. et al, 1999, in preparation (paper II) 

Gerin M., Phillips T.G., Keene J., Betz A.L., Boreiko R.T., ApJ 500, 329 
de Graauw et al., 1996,A&A 315, L49 
Kessler M.F. et al., 1996, A&A 315, L27 

Lemaire J.L., Field D. Gerin M., Leach S., Pineau des Forets G., Rostas F., Rouan D., 
1996, A&A 308, 895 

Martini P., Sellgren K., Mora J.L., 1997, ApJ 484, 296 

Neufeld D.A., Melnick G.J., Harwit M., 1998, ApJ 506, L75 

Salama A. et al., 1997, ESA SP-419, 17 

Skinner S.L., Brown A., Stewart R.T., 1993, ApJSS 87, 217 

Smith M.D., Davis C.J., Lioure A., 1997, A&A 327, 1206 



Sternberg A. & Neufeld D.A., 1999, preprint 

Timmermann R., Bertoldi F., Wright CM., Drapatz S., Draine B.T., Haser L., Sternberg 
A., 1996, A&A 315, L281 

Wesselius P.R. et al, 1996, A&A 315, L197 

Wright CM., Drapatz S., Timmermann R., van der Werf RR, Katterloher R., de Grauuw 
T., 1996, A&A 315, L301 



This manuscript was prepared with the AAS macros v4.0. 



-13- 



Fig. 1.— Spectra of the H2 S(0),S(1),S(2),S(3),S(4) and S(5) towards the peak of the PDR 
associated with NGC 7023 (R.A. (2000): 21^01•^32^5, Dcc(2000): 68°10'27".5). The S(l), 
S(2), S(3), S(4), and S(5) hnes are offset by 4,8,12,18 and 22 W cm'^ /^m"^ respectively. 

Fig. 2. — Rotational diagram of the H2 lines towards peak of the PDR associated with 
NGC 7023. The dashed-line connect the para- H2 levels while the continuous line connect 
the ortho- H2 levels. Note that there is an offset between the two curves and that this offset 
decreases for high temperatures. 

Fig. 3. — Comparison between the observational data (filled squares) and the rotational 
diagram predicted by the model by Burton et al. (1992) for Go = 10"^ and n = 10® cm~^ 
corrected for different values of the OTP ratio, 3 (open triangles), 2 (open circles) and l(open 
squares) . 
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Table 1: Observational parameters 
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